Stability measures in metastable states with Gaussian colored noise 
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We present a study of the escape time from a metastable state of an overdamped Brownian particle, 
in the presence of colored noise generated by Ornstein-Uhlenbeck process. We analyze the role of 
the correlation time on the enhancement of the mean first passage time through a potential barrier 
and on the behavior of the mean growth rate coefficient as a function of the noise intensity. We 
observe the noise enhanced stability effect for all the initial unstable states used, and for all values 
of the correlation time r c investigated. We can distinguish two dynamical regimes characterized 
by weak and strong correlated noise respectively, depending on the value of r c with respect to the 
relaxation time of the system. 

PACS numbers: 05.40.-a, 87.17.Aa, 87.23. Cc, 82.20.-w 
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I. INTRODUCTION 

The problem of the lifetime of a metastable state has 
been addressed in a variety of areas, including first- 
order phase transitions, Josephson junctions, field the- 
ory and chemical kinetics [l], [|J . Recent experimental 
and theoretical results show that long-live metastable 
states are observed in different areas of physics 0,(3. 
Experimental and theoretical investigations have shown 
that the average escape time from metastable states 
in fluctuating potentials presents a nonmonotonic be- 
havior as a function of the noise intensity with the 
presence of a maximum [1, 0, 0j- This is the noise 
enhanced stability (NES) phenomenon: the stability 
of metastable states can be enhanced and the aver- 
age life time of the metastable state increases non- 
monotonically with the noise intensity. This resonance- 
like behavior contradicts the monotonic behavior of the 
Kramers theory The occurrence of the enhance- 
ment of stability of metastable states by the noise has 
been observed in different physical and biological sys- 
tems @, i, i 0, i, M El El El 0, HI . Very recently 
NES effect was observed in an ecological system an 
oscillator chemical system (the Belousov-Zhabotinsky 
reaction) [l7[ and in magnetic systems Interest- 
ingly in Ref. [T3| the stabilization of a metastable state 
due to noise is experimentally detected and a decreas- 
ing behavior of the maximum Lyapunov exponent as a 
function of the noise intensity is observed. 

A generalization of the Lyapunov exponent for 
stochastic systems has been recently defined in Ref. [l9j 
to complement the analysis of the transient dynamics of 
metastable states. This new measure of stability is the 
'mean growth rate coefficient' (MGRC) A and it is eval- 
uated by a similar procedure used for the calculation of 
the Lyapunov exponent in stochastic systems [2fj| . By 
linearizing the Langevin equation of motion (see next 
Eq.[4j, we consider the evolution of the separation Sx(t) 
between two neighboring trajectories of the Brownian 
particle starting at xq and reaching xp 



Sx(t) 



d 2 U(x) 
dx 2 



5x(t) — \i(x,t)8x(t) , 



(1) 



and define Aj(x, t) as an instantaneous growth rate. We 
note that, in Eq. {T]), d 2 U(x)/dx 2 is calculated onto the 
noisy trajectory x [£(£)] [19|]. The growth rate coeffi- 
cient Aj (for the ith noise realization), is then defined 
as the long-time average of the instantaneous Xi coeffi- 
cient over t(xo,x_f) fig. I20I [2l| 



A; 



t(x ,xf) 



■(x ,xf) 



Xi(x, s)ds. 



(2) 



In the limit t(xq,xf) — > oo, Eq. ([2]) coincides formally 
with the definition of the maximum Lyapunov expo- 
nent, and therefore, the Aj coefficient has the meaning 
of a finite-time Lyapunov exponent. This quantity is 
useful to characterize a transient dynamics in nonequi- 
librium dynamical systems [TtI Il9l ] . The mean growth 
rate coefficient A is then defined as the ensemble aver- 
age of the growth rate coefficient Aj 



A = < Ai > 



(3) 
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over the noise realizations. The mean growth rate co- 
efficient has a nonmonotonic behavior as a function of 
the noise intensity for Brownian particles starting from 
unstable initial positions [l9j . This nonmonotonicity 
with a minimum indicates that A can be used as a new 
suitable measure or signature of the NES effect. 

The inclusion of realistic noise sources, with a finite 
correlation time, impacts both the stationary and the 
dynamic features of nonlinear systems. For metastable 
thermal equilibrium systems it has been demonstrated 
that colored thermal noise can substantially modify the 
crossing barrier process Q . A rich and enormous liter- 
ature on escape processes driven by colored noise was 
produced in the 80's [H, [23|, [24|. More recently many 
papers investigated the role of the correlated noise on 
different physical systems [H, [H, HI HI, [H, HI , which 
indicates a renewed interest in the realistic noise source 
effects. 

In this work we present a study of the average de- 
cay time of an overdamped Brownian particle subject 
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to a cubic potential with a metastable state. We focus 
on the role of different unstable initial conditions and 
of colored noise in the average escape time. The effect 
of the correlation time t c on the transient dynamics of 
the escape process is related to the characteristic time 
scale of the system, that is the relaxation time inside 
the metastable state r r . For t c < r r , the dynamical 
regime of the Brownian particle is close to the white 
noise dynamics. For r c > r r , we obtain: (i) a big shift of 
the increase of the average escape times towards higher 
noise intensities; (ii) an enhancement of the value of 
the average escape time maximum with a broadening 
of the NES region in the plane (r, D) , which becomes 
very large for high values of r c ; (iii) the shift of the pe- 
culiar initial position x c (towards lower values), found 
in our previous studies 0, , which separates the set 
of the initial unstable states producing divergency, for 
D tending to zero, from those which give only a non- 
monotonic behavior of the average escape time; (iv) the 
entire qualitative behaviors (i— iii) can be applied to the 
standard deviation of the escape time; (v) the shift of 
the minimum values in the curves of the mean growth 
rate coefficient A; (vi) trend to the disappearance of the 
minimum in the curves of A, with a decreasing mono- 
tonic behavior for increasing r c ; (vii) trend to the disap- 
pearance of the divergent dynamical regime in r, with 
increasing r c . The paper is organized as follows. In 
the next section we introduce the model. In the third 
section we show the results and in the final section we 
draw the conclusions. 



profile has a local stable state at x = and an unstable 
state at x = 1 (see Fig. Q}. The relaxation time for 



the metastable state at x 



d z U(x) 
dx 1 



2a, 



x=0 



is T r 

which is the characteristic time scale of our system. For 
our potential profile we have r r = 0.6. 
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FIG. 1: (Color online) The cubic potential U(x) = ax 2 — bx :i 
with the various initial positions investigated (dots) , namely 
x = 1.2, 1.3, 1.4, 1.5, 1.6. The parameters of the potential 
are: a = 0.3,6 = 0.2. For the white noise case, x c — 1.5 
is the critical initial position which separates the set of the 
initial unstable states producing divergency, for D tending 
to zero, from those which give only a nonmonotonic behavior 
of the average escape time 



ve only ; 

s ta- 



il. THE MODEL 

The starting point of our study is the Langevin equa- 
tion 

* = -^M + , W (4) 
where r\ (t) is the Ornstein-Uhlenbeck process 

drj = -h]dt + kVD dW(t) (5) 

and W(t) is the Wiener process with the usual sta- 
tistical properties: (£(£)) = and (£(i)£(i + t)) = 
8{t). The system of Eqs. dU) and ([5]) represents a two- 
dimensional Markovian process, which is equivalent to a 
non-Markovian Langevin equation driven with additive 
Gaussian correlated noise, with ry(t) obeying the follow- 
ing statistical properties (??(£)) = and (r](t)rj(t + t)) = 
(kD/2)e~ kT , for t — ► oo and 77(0) = 0. Here l/k = r c is 
the correlation time of the process. The integration of 
Eq. ([5]) yields in the limit r c — > the white noise term 

ft e -(t-t')/r c 

lim r](t) = 2VD / lim dW(t') = y/Dtft), 

t c ^o J t c ^o 2r c 

and the stationary correlation function of the Ornstein- 
Uhlenbeck process gives in the limit r c — > the cor- 
relation function of the white noise: lim Tc ^o(ri(t)'r}(t + 
t)) = DS(t). The potential U(x) used in Eq. @ is 
U{x) = ax 2 - bx 3 , with a = 0.3, b = 0.2. The potential 



III. RESULTS 

The calculations of the average escape time as a func- 
tion of the colored noise intensity have been performed 
by averaging over TV = 20, 000 realizations the numer- 
ical solution of the stochastic differential equation ([4]). 
The absorbing boundary for the escape process is put 
on if = 20, and the maximum simulation time is 
T m ax = 10,000 a.u.. For all the initial unstable states 
(see Fig. [IJ and all the correlation times considered we 
find an enhancement of the mean first passage time 
(MFPT) t with respect to the deterministic time. 

In Fig. [2] the calculations performed with low colored 
noise (r c = 0.01) for the mean first passage time r and 
the mean growth rate coefficient A are shown. We see as 
signatures of the NES effect a maximum in the curve of 
t and a minimum in that of A. In the inset of Fig. [2^ the 
standard deviation of the first passage time as a func- 
tion of noise intensity is reported. We note that the 
behaviors of r and A in this low colored noise regime 
(t c = 0.01) coincides with those obtained in the white 
noise case 19]. Moreover by comparing the theoretical 
predictions of r (see Eq. (3) of Ref. [lj|) with direct 
numerical simulations of the Langevin equation, a very 
good agreement is obtained (see Fig. 3 of Ref. [lj|). 
In Fig. [3] the semi-Log plots of the fraction of particles 
Ni/N reaching the threshold position x t = 0.5 into the 
potential well, within the T max , as a function of noise 
intensity D, with the same initial conditions of Fig. [1] 
are shown. This threshold position xt corresponds to 
the concavity change of the potential and is considered 
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FIG. 2: (Color online) Panel a: Log-Log plot of the mean 
first passage time r as a function of noise intensity D in 
the case of correlated noise with r c = 0.01, for the four 
initial positions investigated (see Fig. [TJ. Inset: the related 
standard deviation as a function of the noise intensity D. 
The dotted straight line at D — D 3 separates the simulation 
data representing the Brownian particles escaped within the 
maximum simulation time T max for D > D a , from those 
representing the particles partially trapped within the well 
for a time greater or equal to T max for D < D s . Panel b: 
Mean growth rate coefficient A as a function of the noise 
intensity D, with the same initial positions of Fig. [1] 




FIG. 3: (Color online) Semi-Log plot of the fraction of parti- 
cles Ni/N reaching the threshold position x t = 0.5 into the 
potential well, within the Tmax, as a function of noise inten- 
sity D. This threshold position xt corresponds to the flex 
point of the potential, where the instantaneous growth rate 
Xi(x,t) is equal to zero. The correlation time of the noise is 
t c — 0.01., with the same initial conditions of Fig. [T] 



for this reason as a reference indicator for the effective 
entrance of the particle into the well. It is possible to 
observe that for very low noise intensity none particle 
enters into the well within the T max considered, and the 
estimation of the stability measures take their deter- 
ministic values. We note that the behavior of the mean 
growth rate coefficient as a function of the noise inten- 
sity is strongly affected by the characteristic potential 
shape of a metastable state. The curves shown in Fig. [3] 
clarify the behavior of A in the limit of D — > 0. In fact 
the position xt — 0.5 is the flex point of the potential, 
where the instantaneous growth rate \{x, t) is equal to 
zero. We see that for low noise intensities the fraction 
Ni/N goes to zero, producing an increasing behavior of 
the MGRC (see Fig.rjb). 

The behaviors of the MFPTs as a function of the 
noise intensity D with other values of r c are shown in 
Fig. 2J We clearly observe two dynamical regimes de- 
pending on the value of r c with respect to the relaxation 
time of the system (r r = 0.6): (a) weak colored noise 
(0 < t c < r r ) and (b) strong colored noise (r c > r r ). By 
observing Fig. 0Ji (i~ c = 0.1) we can see that the quali- 
tative behavior of MFPT shown in Fig. [2^ is recovered. 
In the weak color noise regime we can still observe the 
divergent behavior of MFPTs for x max < xo < x c and a 
non monotonic behavior for xq > x c , with x c = 1.5. By 
increasing the value of the correlation time (r c > r r ) we 
observe a large displacement of the maximum of MFPT 
towards higher values of noise intensity and a shift of 
the peculiar initial position x c towards lower values. For 
t c = r r = 0.6, x* ~ 1.4, and for r c = 1, x* ~ 1.3 (see 
Figs. [4)d and HJ:), where x* c is the peculiar initial posi- 
tion of the Brownian particle in the presence of colored 
noise. We note that x c = 1.5 is a fixed value for white 
noise case [l9| , while the position x* is a variable quan- 
tity for colored noise and it is depending on the value of 
the correlation time of the noise. Moreover, we observe 
a broadening of the NES region, which becomes very 
large for high values of the correlation time r c . The 
NES region is the area where enhanced stability of a 
metastable state is observed. In other words it is the 
area under each curve of r vs D (see Figs. [5^ and |3|, 
where the values of r are greater than the deterministic 
dynamical time related to the particular initial position 
investigated (see also Fig. 1 in Mantegna and Spagnolo, 
1998, Ref. [|). 

The asymmetry of the potential profile with respect 
to the x coordinate makes more effective the correla- 
tion of the noise for Brownian particles moving from 
left to right. This means that, at very low noise inten- 
sities of the colored noise, the particles inside the po- 
tential well will escape more easily with respect to the 
white noise case. Therefore, the trapping effect, which 
is responsible for the divergent behavior for any initial 
unstable state within the range x max < x D < x c will 
happen in a restricted range of initial positions, that is 
Xmax < x a < x* c with x* c < x c . Specifically this pe- 
culiar position x* is shifted towards decreasing values 
of the x coordinate for increasing correlation time r c of 
the noise source. In Fig. [2^, and all panels of Fig. [4] the 
dotted straight line at D = D s separates the simula- 
tion data representing the Brownian particles escaped 
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FIG. 4: (Color online) Log-Log plot of the MFPT r as a function of noise intensity D for the same initial positions of 
Fig. [1] and for different values of the correlation times t c , namely r c = 0.1,0.6, 1, corresponding respectively to the weak, 
intermediate and strong colored noise dynamical regime. The dotted straight line at D — D a separates the noise values for 
which all the Brownian particles escape (D > D 3 ), from those for which the particles are partially trapped into the well 
within the T max (D < D s ). 




FIG. 5: (Color online) Log-Log plot of the standard devia- 
tion o as a function of the noise intensity D for r c = 0.6 and 
the same initial positions of Fig. [1] The dotted straight line 
indicates the value of the noise intensity D s which separates 
the simulation data representing the trapped Brownian par- 
ticles from those escaped within T ma x- Inset: Log- Log plot 
of the ratio cr/r as a function of noise intensity D. 



within the maximum simulation time T max for D > D s , 
from those representing the particles partially trapped 
within the well for a time greater or equal to T max for 
D < D s . This means that the simulation data obtained 
for D < D s underestimate the real data in the divergent 
dynamical regime. In fact if we prolong the maximum 
simulation time T max we obtain more approximate val- 
ues for r and a and the divergent behavior will be visible 
at lower noise intensities. As a consequence D s will be 
shifted towards lower values. 

For high values of the noise intensity all the plots 
show a monotonic decrease behavior as a function of 
noise intensity collapsing in a unique curve. Moreover 
the slope of this limit curve becomes flatter by increas- 
ing the correlation time. This means that the NES ef- 
fect involves more and more orders of magnitude of the 
noise intensity. The effect of the colored noise is there- 
fore to delay the escape process or in other words to 
enhance more and more the stability of the metastable 
state for increasing values of the noise intensity. 

In Fig. [5] the standard deviation a of the first pas- 



sage time distribution for r c = 0.6 is shown. We see a 
huge increase of the a for low values of noise intensity, 
demonstrating a strong enlargement of the distribution 
when the particle feels a noise intensity comparable with 
the height of potential barrier. Similarly to the MFPTs, 
color induces a shift in the divergent behavior of a. The 
relative measure of the width with respect to the mean 
value is shown in the inset of Fig. [5J where the ratio 
cr/r is plotted. This ratio reveals a nonmonotonic be- 
havior with a minimum, demonstrating the existence of 
a noise intensity for which the width of the first passage 
time distribution is the minimum related to its mean. 
In other words this value corresponds to a maximum of 
precision in the measure of r. This optimal noise in- 
tensity is shifted toward high noise values by increasing 

The behavior of the mean growth rate coefficient A 
as a function of the noise intensity D for different val- 
ues of the noise correlation time is shown in Fig. [6l In 
the weak color noise regime we observe a nonmonotonic 
behavior with a minimum for all the initial positions 
investigated with a shift in the position of the mini- 
mum towards higher noise intensities. In the strong 
color regime the minimum, which represents a trapping 
phenomenon for a finite time, is visible for the diver- 
gent behavior of MFPTs for x max < xq < x* c and it 
is shifted towards higher noise intensities by increas- 
ing the correlation time. For initial positions xq > x* c , 
the minimum tends to disappear, but at the same time 
the A parameter decreases monotonically with increas- 
ing noise intensity, showing a trapping phenomenon at 
higher noise intensities. This trend to the disappear- 
ance of the minimum in the curves of A, corresponds to 
the restricted range of the initial positions for which we 
observe a divergent behavior of r, that is to the trend of 
disappearance of this divergent behavior. We note that 
the behaviour of A as a function of the noise intensity 
D obtained in our analysis is in qualitative agreement 
with that obtained by the experimental investigation 
of the stabilization of a metastable state in an oscilla- 
tory chemical system (the Belousov-Zhabotinsky reac- 
tion) [17[. Specifically the decreasing behavior of the 
maximum Lyapunov exponent of Fig. 2 of Ref. [17l | is 
in qualitative good agreement with the behavior of the 
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FIG. 6: (Color online) Semi-Log plot of mean growth rate coefficient A as a function of noise intensity D for the same initial 
positions x and the same values of the correlation times r c of Fig. [4] namely r c = 0.1, 0.6, 1, corresponding respectively to 
the weak, intermediate and strong colored noise dynamical regime. 




D 

FIG. 7: (Color online) Semi-Log plot of the fraction of par- 
ticles Ni/N reaching the threshold position x t = 0.5 into 
the potential well, within the T ma x, as a function of noise 
intensity D. This threshold position xt corresponds to the 
flex point of the potential, where the instantaneous growth 
rate \%(x,t) is equal to zero. The correlation time of the 
noise is r c = 0.6, with the same initial conditions of Fig. [T] 

curves in Fig. [SJa and Fig. This could be ascribed 
to the correlation time always present in noise sources 
used in any experimental set-up. 

In Fig. [7] we report, for all the initial positions in- 
vestigated and for r c = 0.6, the semi-Log plot of the 
fraction of particles Ni/N entering into the potential 
well up to the position x t = 0.5, within the T max , as a 
function of noise intensity D. At very low noise intensi- 
ties and for increasing values of the correlation time r c , 
the particles have difficulty to enter into the potential 
well, within the T max considered, shifting the entrance 
statistics towards higher values of the noise intensity. 



IV. CONCLUSIONS 



In this work we analyzed the effect of the colored 
noise, generated by an Ornstein-Uhlenbeck process, on 
the enhancement of the mean first passage time in a 
cubic potential with a metastable state and on the min- 
imum of the mean growth rate coefficient as a function 
of the noise intensity. We analyze different initial un- 
stable states. We obtain NES effect for all the initial 
positions investigated and an enhancement of the NES 
region for increasing values of correlation times. The re- 
sults obtained for a particle moving in a cubic potential 
are quite general, because we always obtain NES effect 
when a particle is initially located just to the right of 
a local potential maximum and next to a metastable 
state, in the escape region. 

In experiments real noise sources are correlated with a 
finite correlation time. As a consequence the NES effect 
can be observed at higher noise intensities with respect 
to the idealized white noise case. The enhancement and 
the shift of the NES region, towards higher values of the 
noise intensity, allows to reveal experimentally the NES 
effect only by using a suitable correlation time r c in the 
noise source. 

This work was supported by MIUR. 
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